Numerical study of domain coarsening in anisotropic stripe patterns 
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We study the coarsening of two-dimensional smectic polycrystals characterized by grains of obhque 
stripes with only two possible orientations. For this purpose, an anisotropic Swift-Hohenberg equa- 
tion is solved. For quenches close enough to the onset of stripe formation, the average domain 
size increases with time as t^^'^ . Further from onset, anisotropic pinning forces similar to Peierls 
stresses in solid crystals slow down defects, and growth becomes anisotropic. In a wide range of 
quench depths, dislocation arrays remain mobile and dislocation density roughly decays as t~^^^ , 
while chevron boundaries are totally pinned. We discuss some agreements and disagreements found 
with recent experimental results on the coarsening of anisotropic electroconvection patterns. 

PACS numbers: 47.54.-|-r, 61.30.Jf, 64.60.Cn 
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I. INTRODUCTION 

Coarsening occurs when a system is rapidly quenched 
below a transition point into a phase with broken symme- 
tries. The spatio-temporal evolution following a quench 
is relatively well known when the broken symmetry phase 
is characterized by a local order parameter that is spa- 
tially uniform (e.o. the local magnetization of a ferro- 
magnetic phase) Q, . The growth of spatial correla- 
tions, driven by domain growth or the annihilation of 
topological defects, usually obeys a dynamical scaling re- 
lation and the correlation length, or "domain size" R, 
grows as a power law of time with a well defined expo- 
nent. Classification schemes have been established for 
the main different cases |3, Q ■ 

The situation is much less understood for phases char- 
acterized by a local order parameter that is spatially 
modulated. Systems forming periodic patterns (stripes, 
hexagons) with a well defined periodicity can be ob- 
served in numerous physical systems, such as Rayleigh- 
Benard convection, diblock-copolymer melts, magnetic 
materials, or Turing reaction-diffusion systems 0|. After 
a quench into a stripe phase, two dimensional config- 
urations are composed of many domains differently ori- 
ented, including grain boundaries, dislocations and discli- 
nations. Numerical I|JiJlJl |1 El El IH El as well 
as experimental Il5l Il6[ll7l | studies have established 
that it is difficult, if not impossible, to reduce the or- 
dering dynamics of stripes to one of the class known for 
uniform phases. There is still some debate regarding the 
growth mechanisms, the value of the growth exponent, 
whether dynamical scaling holds or not, or whether the 
system may involve various characteristic length scales 
growing with different exponents. In contrast with uni- 
form phases, the coarsening rates depend significantly on 
the quench depth. Far away from the bifurcation thresh- 
old of stripe formation (large quenches), numerical solu- 
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tions of the Swift-Hohenberg equation show that coarsen- 
ing stops at large time, i.e. the system remain frozen in 
macroscopically disordered configurations H • Based 
on an analysis of the law of motion of a grain bound- 
ary through curved stripes, it was recently proposed that 
a single growth exponent could be introduced, but for 
vanishingly small quenches only IT^ . From dimen- 
sional arguments, a i? ~ t^^^ growth law was derived in 
that regime, in good agreement with numerical results 
at small quenches jO; E^ ■ The freezing observed at 
finite quenches was attributed to the presence of a peri- 
odic pinning potential (generated by the pattern itself) 
acting on grain boundaries. 

In the present paper, we consider a closely related 
problem where similar questions remain open, and that 
has not been investigated numerically so far: the coarsen- 
ing of anisotropic stripe patterns. Oblique rolls making 
only two possible angles {9 or —6, fixed) with respect to 
a particular axis can be observed in electroconvection of 
nematic liquid crystals . Studying the ordering dynam- 
ics on this system is motivated by various reasons. First, 
one can intend a comparison with available experimental 
data, since coarsening experiments have been recently 
conducted in electroconvection 0, Second, the 

polycrystalline structures of oblique stripes have a rel- 
atively simpler geometry than those of isotropic stripes: 
the constraint of the fixed angle prevents the formation 
of disclinations. Therefore, the topological defects are es- 
sentially dislocations and (chevron) grain boundaries sep- 
arating domains differently oriented. This situation can 
be seen as a smectic analogue of the structures formed 
by grains in polycrystalline solids, where disclinations are 
also absent 18]. 

We consider in the following an extension of the Swift- 
Hohenberg equation for oblique stripes in two spatial di- 
mensions. This model, proposed by Pesch and Kramer 
for describing electroconvection E3 , is recalled in Section 
IhI In Section [nil we investigate small quenches: the nu- 
merical results show that coarsening is driven by surface 
tension and a growth law t^^"^ is observed for various char- 
acteristic length scales, like in Model A Q. The results 
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qualitatively change at larger quenches fSection llV|l : the 
characteristic length scales associated with dislocations 
and chevron boundaries start to evolve differently, and 
the associated effective growth exponents progressively 
decrease as the quench depth increases. However, the ef- 
fective exponent of the dislocation density remains fairly 
constant for a relatively wide range of quench depths. 
This feature can be explained by the fact that disloca- 
tions have a much lower pinning potential than chevron 
boundaries. We qualitatively justify this feature from 
weakly nonlinear analysis arguments. The dislocation 
exponent is close to the value of 1 /3 in that intermediate 
range, in agreement with the value measured in recent 
experiments . Some conclusions are presented in Sec- 
tion El 



II. MODEL EQUATION 

Electroconvection in nematic liquid crystals is a 
paradigm of anisotropic pattern formation |2(TI |. If a 
nematics is placed between two glass plates properly 
treated, its director can be aligned along a preferential 
direction, say the a;-axis. When an external a.c. elec- 
tric field is applied in the direction normal to the plane, 
periodic rolls appear above a threshold. As the voltage 
is increased (the frequency being fixed in some proper 
range), bifurcations to various phases can be observed: 
"normal" rolls, with a wave vector directed along the x- 
axis, usually appear first. This phase can be followed by 
a transition to "oblique" rolls, of interest here, character- 
ized by a wavevector with two possible orientations with 
respect to the x-axis, and ~9. 

Although the theoretical understanding of electrocon- 
vection patterns based on constitutive equations is still 
incomplete, some nonlinear models that rely on equations 
for a local order parameter and on symmetry arguments 
have been proposed. Some time ago , Pesch and Kramer 
introduced an anisotropic model [l^] that exhibits a tran- 
sition from normal to oblique rolls: 
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d/dx{y). In equation ip{x,t) is a lo- 



dt 
with 

cal dimensionless order parameter, interpreted as a small 
lateral elastic displacement; c and rj are dimensionless 
constants modeling the loading forces and anisotropic 
bending constants; kg is the wavenumber of the base pe- 
riodic pattern, and C a "coherence" length that will be 
set to 1/ko for simplicity here. The dimensionless pa- 
rameter r is chosen as the main control parameter. For 
c = 77 = the above equation reduces to the well known 
Swift-Hohenberg model of Rayleigh-Benard convection (r 
is the reduced Rayleigh number in that case) . The model 
derives from a Liapunov "free-energy" functional. It 
can be recast as 



with 

-2r,{d,dy^f + c{dl^f] . (3) 

The functional F monotonically decreases with time, 
dF/dt < 0. 

Linear stability analysis of Eq.QJ around the state 
Tp^x^t) = shows that two modes of finite wavenum- 
ber k ~ px + qy {x and y are unitary vectors) become 
marginally unstable when the control parameter r in- 
creases and crosses some threshold values and re"-* : 
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-ic - .0 ~^ (obhque rolls) (4) 
c + 277 — 7;'^ 

r(") =0, {pl^kl, ql = Q} (normal rolls) (5) 

In the above relations, we have considered the case c > 0, 
the condition for which the instability at Tc = is to 
normal rolls (pc = 0). In the oblique phase, rolls make 
an angle 6 — ± arctan[-\/77/ (c + 77)] with the y-axis. As 
obvious from Eq.l^J, oblique rolls can only be observed 
for 77 > 0. Hence, when rj is tuned from negative to 
positive values, a transition form normal to oblique rolls 
can occur. In order to study the coarsening of oblique 
rolls, we will chose c > and 77 > in the following. 

Oblique rolls have the lowest threshold value (rc°'' < 0) 
and we re-note Vc = for simplicity. As already noted 
by Pesch and Kramer, the structures that are likely to 
be observed for r in the range [re, 0] may not be oblique 
rolls, but more complicated, nonlinear structures ("un- 
dulated" rolls) that are not of interest here. When nu- 
merically solving Eq.(^) with random initial conditions 
and < r < 0, we actually observed that oblique rolls 
never appeared. On the other hand, configurations of 
oblique rolls very similar to those observed experimen- 
tally ^1 are always observed when setting r > 
instead. Normal rolls were never observed in the runs 
presented in the following Sections. This is consistent 
with the weakly nonlinear analysis of Eq.lQJ that pre- 
dicts that oblique rolls have a lower free-energy F [l9l| . 
Therefore, we define the quench depth as 

e = r - Tc, (6) 

with Tc given by Eq.Q). We always chose e larger than 
\rc\ {r > 0). 



III. COARSENING KINETICS NEAR ONSET 

« 1) 

We numerically solve Eq. JQ) by using a pseudo-spectral 
method and a time integration procedure whose descrip- 
tions can be found in Ref.(21l]. The space is discretized 
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FIG. 2: Defects (marked in black) of a configuration with 
same parameters as in Fig0 at a larger scale. 

on a square lattice of 1024^ nodes with a lattice size Ax 
set to unity. The base period Aq = 27r/fcQ of the pattern 
is fixed to 8Aa;. The time integration scheme is stable for 
relatively large value of the time step, which is fixed to 
0.5 in dimensionless time units. The initial condition for 
■f/i is a random field with Gaussian distribution, of zero 
mean and variance -y/e/S. 

Figure ^ displays in gray scale the order parameter ip 
at time t = 1200 time units, for a run with c = 12 and 
ry = 0.5 (the angle of the rolls with the vertical axis is 
11.31°). The quench depth is small, and has been set 
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FIG. 3: Defect density as a function of time. From bottom 
to top. (+) symbols : c = 12, 77 = 0.5 {6 = 11.31°), e = 
1.9|rc| = 0.0372. (0) symbols : c = 6, 77 = 0.25 (6* = 11.31°), 
e = 1.9|rc| = 0.0184. (Q) symbols: c = 3, »? = 0.25 {6 = 
15.50°) , e = 1.9|rc| = 0.0345. Average are performed over 13 
runs in each case. Solid lines are guides to the eye. 



to e = 1.9|rc| ~ 0.0372. The configuration is that of a 
smectic polycrystal: Most of the defects present are grain 
boundaries separating zig and zag rolls, and few isolated 
dislocations can be observed. Due to the asymmetry of 
the problem, two kind of boundaries can be roughly dis- 
tinguished: "horizontal" chevron boundaries where the 
roll orientation changes rather smoothly from one grain 
to the other, and "vertical" (or inclined) boundaries, that 
are made of dense arrays of dislocations. This distinction 
is not very sharp, as one go continuously from one situ- 
ation to the other, corresponding to boundaries of "low" 
and "high" dislocation density, respectively. The defect 
field shown in Figure |21 is obtained from FigQ] by using 
a Fourier filtering procedure. The dark areas correspond 
to defected regions. Dislocations tend to be distributed 
along string-like structures, like in crystals. 

At very large times, grain boundaries are weakly 
curved and isolated dislocations lying inside a grain are 
relatively rare. These observations agree with recent elec- 
troconvection experiments [l5| |. where a mechanism for 
the formation of isolated dislocations was identified: A 
shrinking bubble can be roughly pictured as delimited 
by two vertical and two horizontal grain boundaries. In 
some cases, the two vertical boundaries are not composed 
by the same number of dislocations, therefore, some dislo- 
cations can not annihilate with others of opposite Burgers 
vector when the bubble shrinks. However, this situation 
occurs rarely. 

We study the time evolution of the defect density pd, 
defined as the fraction of area occupied by the black re- 
gions in Fig. |2] We perform three series of runs at small 
quenches, each satisfying e — l.Olrd for different choices 
of the parameters {c,?7,r} (e = 0.0372,0.0184,0.0345 
respectively). Figure 13 shows a summary of the data ob- 
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FIG. 4: Relative Liapunov free-energy per unit area (^4 is the 
system area) as a function of time. The legend is the same as 
in Figure 121 Solid lines are guides to the eye. 

tained. In each cases, the results are consistent with the 
law 

P,-t-'/', (7) 

which corresponds to a defect characteristic length scale 
growing as t^/^. This result seems to be fairly indepen- 
dent of the angle 6. 

We next investigate the time evolution of the Liapunov 
functional F given by Eq.Q. If Fq denotes the value 
of F for a perfectly ordered system, then the quantity 
AF = F — Fq represents the excess energy due to defects. 
The system free energy F decreases like the total length 
of grain boundaries. From Q, one should expect 

F-Fo-t"^/2. (8) 

The numerical data plotted in Figure 0] for small values 
of e are consistent with this scaling relation as well. 

The time evolution of many coarsening systems is 
self-similar: the large scale structure of successive con- 
figurations is statistically time invariant, provided that 
spatial variables are rescaled by a proper length. We 
thereby analyze the structure factor, defined as the 
Fourier transform of the two-point correlation function, 
S{k) = {ip{k, t)ip{—k, t)), the braces representing average 
over initial conditions. S{k) is maximum for any of the 
(four) wavevectors kz characterizing zig and zag rolls. At 
any given time, we numerically observe that 5' is maxi- 
mal for kz ~ kc, given by Eq.l^J. Therefore, the selected 
stripe periodicity and orientation in the polycrystalline 
structure are that of the marginal wavevectors kc deter- 
mined from linear stability analysis. A similar situation 
is encountered for the isotropic Swift-Hohenberg model, 
and is thought to be more generally a property of poten- 
tial systems (7|. Near any peak kz of S, we propose the 
following scaling ansatz 

5(Jfc||,5fci,i)=e||We±(i) (9) 
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FIG. 5: Structure factor (5fc|| , t) at four different times: 
t = 510^(0), 210^(*), 510'^{+), 210^(A). At each time, the 
curve has been rescaled according to Eq.l^JJ, where £,\\{t) is 
defined as S{Sk^^ — 0,t). The parameters are c = 6, = 
0.25, e — 1.9\rc\. (Averages over 40 independent runs.) 
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FIG. 6: Length determined from the maximal intensity of 
the structure factor. See Figure |21 for legends. The solid line 
is a guide to the eye. 

where Sk = k — kz = + Sk±k±, with fc|| and 

k±_ denoting the unit vectors longitudinal and transverse 
to the wavevector kz, respectively. f{x) is a scaling 
function, S.\\(t) and ^±{t) are a priori two characteris- 
tic lengths describing grain growth in the directions nor- 
mal and parallel to the rolls, respectively. Let us define 
£"11 (J/eii , t) = Sk±S{Sk,t), a symmetric relationship 
defining S±. From 0, one obtains the scaling ansatz 

£||(Jfc||,t) = e||(0 5[e||(0'5fc||] • (10) 

Figure displays S^^ as a function of (5fc|| at various 
times (t = 510^, 210^, 5 10^ 510''), for {c = 6, = 
0.25, e = 1.9|rc|}. The scaling relation H10|) holds over 
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nearly two decades, despite of a slight widening at large 
times. The length £,\\(t), taken from the maximum value 
of S*!! is plotted as a function of time in FigEl The results 
are in good agreement with 

m-t'^'- (11) 

We have not found ^±{t), determined from 5"^, to be a 
convenient length scale to characterize coarsening. ^± is 
larger than ^|| by a factor varying between 5 (at short 
times) and 2 (at large times). This could be due to some 
phase correlations of longer range than those associated 
with domain walls; these features were not investigated. 

The above results are consistent with a coarsening pro- 
cess driven by grain boundary surface tension and in- 
volving a single characteristic length scale, like for the 
dynamics of Model A for a non-conserved order parame- 
ter Q. The law R ~ t^/^ is also expected to describe the 
kinetics of grain growth in solid polycrystals [22|. This 
situation differs markedly from isotropic stripes (where 
R ^ t^/^ milUlll)- Note that the available experi- 
mental studies on anisotropic stripes have reported much 
slower coarsening laws than Eq. namely t^^^ or t^/^ 

EHIl. 

IV. COARSENING AT LARGER QUENCH 
DEPTHS 

A. Phenomenology of pinning forces in modulated 
phases 

In solid crystals, for an isolated dislocation to glide 
from one raw of atoms to the next one, there is a fi- 
nite energy cost corresponding to the raws that have to 
be compressed or dilated during the move. The result- 
ing elastic force, the Peierls stress, tends to prevent the 
glide of dislocations ^3^2^- Therefore, dislocations tend 
to be pinned in positions (periodically spaced) of mini- 
mum energy, and motion takes place only if an external 
stress larger than the critical Peierls stress is applied. Re- 
markably, the defects present in systems that form peri- 
odic patterns are also subjected to similar pinning forces. 
Their origin is nonlinear in that case, and is due to the 
apparition of "non-adiabatic" terms in weakly nonlinear 
expansions. Studies on the Swift-Hohenberg model have 
shown that the laws of motion of grain boundaries (or 
dislocation arrays) involve short range, spatially periodic 
pinning forces 12, 25, 26, 27] . Either in crystals or in 
patterns, the law of motion of a defect takes the general 
form 

li~^v = fi~^dx/dt = f — pcos{kpX), (12) 

where v is the defect velocity, x its position (for a grain 
boundary, the coordinate normal to the interface), fi a 
mobility, / an external force per unit length {e.g, the 
driving force for coarsening) , and p the magnitude of the 
pinning force, that oscillates with a periodicity 27r/fcp 



proportional to the crystal (or base pattern) periodicity 
Ac. Peierls-like pinning forces are usually much smaller 
than the other characteristic elastic forces (like the criti- 
cal threshold shear stress fcr), and have the same general 
approximate expression, valid both for solids |i23i.i24j and 
nonlinear patterns [l2, 12^ : 

p/U^eM-aW/Xc], (13) 

where W is the width of the defect (see further Fig. [SJ, 
and a a constant of order unity. For instance, grain 
boundaries separating domains of stripes have a width 
W ~ Ac/Ve |23|. Therefore, close to the onset of the su- 
percritical bifurcation (e 0^), W becomes very large 
and the pinning potential H13|) can be neglected. On the 
other hand, as the quench depth e increases, W decreases 
and pinning forces can become large enough to affect 
qualitatively defect dynamics. 

During a coarsening processes driven by surface ten- 
sion, the average force / in Ea. H12l) is time dependent: 
/ ~ 7K ^ "f/R(t), with 7 and k the typical interface sur- 
face tension and curvature, respectively. At short times, 
domains have small sizes R, and / is large compared 
with p in Eq. H12|) . Defects move easily and the aver- 
age grain size grows. As a result, the driving force / 
decreases with time. At some point, / may eventually 
become lower than the typical value of p (which is time 
independent). In this case, boundaries become pinned at 
one of the discrete stable positions Xp such that v — Q 
in Ea. (|12|l . This situation is easy to observe numerically 
for large quenches {p large), where asymptotic patterns 
remain only partially ordered ( "glassy" ) , with many im- 
mobile defects 0, [3| ■ On the other hand, in the regime 
e — 0, defect pinning is negligible during the numeri- 
cal time scales studied, and coarsening dynamics is more 
likely to be self-similar and characterized by well defined 
exponents. 

For intermediate quenches, well before all grain bound- 
aries of a system become pinned, the Peierls-like barriers 
are believed to slow down the ordering kinetics. At inter- 
mediate times, one may still be able to fit in some cases 
the numerical results with power laws, R ~ t^^^ . z* now 
represents an effective growth exponent, and z* > z\e^o. 

B. Pinning of anisotropic stripes 

To check whether the arguments presented above ap- 
ply to anisotropic patterns as well, we have performed 
calculations for deeper quenches than in Sec lIIII 

Figure [7^ shows in gray scale the field tp oi a frac- 
tion of the system, obtained for a "moderate" quench 
of depth e = ll|rc| = 0.215, at time t = 500 with 
{c=12,r7 = 0.5}. The two classes of defects previously 
mentioned, the horizontal chevron grain boundaries and 
the dislocations, can now be clearly distinguished. (The 
dislocation field is shown in Fig[7|D.) The chevron bound- 
aries are fairly straight. They remain practically immo- 
bile during the whole coarsening process, that is driven by 




FIG. 7: a) Local order parameter in gray scale, c = 12, 77 = 
0.5, t = 500. The quench is "moderate": e = ll|rc| = 0.215. 
The chevron boundaries are now straight and pinned. 6) Dis- 
location field of a) (same scale) obtained with the Fourier 
filtering procedure. 

dislocation motion onh^This feature was observed in ex- 
periments as well 0, 01 . A detail of a large time config- 
uration {t = 50000) is shown in figure |S1 Like for shallow 
quenches, dislocations tend to organize along string-like 
structures that are generally curved. We interpret the 
immobility of the chevron boundary as caused by strong 
pinning forces. Given two domains of zig and zag rolls, 
the stable positions of a chevron boundary are imposed 
by the phase of the local order parameter, which does not 
change across the boundary when one follows a given roll. 
On the other hand, dislocations are much more mobile, 
suggesting that their pinning potential p is very low, and 
therefore strongly anisotropic with respect to the grain 
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3. 8: Same parameters as in Fig. at f = 50000 (detc 
3 chevron boundaries and dislocations have a width 
1 Wdis respectively. 
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FIG. 9: Relative free energy per unit area as a function of 
time for c = 12, 77 = 0.5 and various quench depths. From 
bottom to top: e = 0.037; 0.098; 0.294; 0.607; 1. 

boundary orientation. 

In the following, we define the effective exponents as- 
sociated to dislocations, free-energy and structure factor, 
respectively: 

Pd^s - t-^'^'^", AF ~ t-'/''-, Cwit) - t^^^'K (14) 

where pdis is the dislocation density, and is determined 
the same way as pd in Section IIIII (the fraction of black 
area of Fig[7|D) 28] . The other quantities have been de- 
fined in Section Inil 

We have plotted the time evolution of the relative free- 
energy per unit area AF in Fig|3 for different quench 
depths. Similar curves are obtained for the dislocation 
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density pjjjg fFig llU|) . and the correlation length ^ (not 
shown). Provided that e < 35|rc|, the curves are still 
reasonably well fitted by power-laws during the first few 
decades considered in the numerical calculations. At 
larger quenches, they rapidly saturates a finite values, 
indicating defect pinning. In these cases, we define the 
effective exponent (arbitrarily) as given at time t — 1000: 



d\n Pdi. 



d\nt 



t=1000 



(15) 



(and similar relations for z*p and z*g.) 

As expected from the discussion of Section FlV Al the 
ordering kinetics slow down noticeably as e is increased. 
All exponents ^^j^, z*p and z*g increase with e. Figure 
111! displays the variations of the different effective expo- 
nents as a function of e. The dislocation exponent z^^g 
differs noticeably from z*p (and z*g): z^^^ > Zp > z*g. 
The behavior of z*p is characterized by two regimes: At 
moderate quenches, Zp gradually departs from z*p — 2 
and slowly increases with e up to a value close to 3. For 
e > 35|rc|, z*p then increases more sharply, the signature 
of a sudden increase of pinning effects. A similar behav- 
ior (although less pronounced) is observed for z*g. The 
behavior of 2^^^ with the quench depth is more abrupt. 
The variations of z^^^ are quite important for small and 
large values of e. The most striking feature is the pres- 
ence of a fairly long plateau at intermediate quenches 
(5|rc| < e < 35|rc|) where z*^^^ remains practically con- 
stant, z^jj, ~ 3. This result agrees with the experimental 
results of 15], where a law t~'^/'^ was reported for the 
dislocation density. 



C. Discussion 

We sketch a possible interpretation of part of the above 
observations, based on the fact that pinning effects are 



strongly anisotropic. We saw that, away from onset, 
chevron boundaries become totally pinned. Meanwhile, 
dislocations are mobile and may still have a very low pin- 
ning potential. From the standard relation for the width 
of a grain boundary in stripe patterns, derived from 
weakly nonlinear analysis [23, let us assume that the 
width of chevron boundaries (Wc) and that of a roughly 
vertical dislocation array iWdisi see Fig. |Sl are given by 
Wc,dis = 5c,dis A/v^, with 5c and 6dis two constants of 
order unity. From observations, let us assume that dislo- 
cations have a larger width than chevrons: Wdis > Wc, 
i.e. Sdis > 5c. 

From the general relation IjlSI) . the pinning potentials 
have an non-analytical behavior at low e, and saturate 
at large t {p ^ exp[— cst/y^). If 5dis > <5c, the pin- 
ning potentials Pc and pdis for chevrons and dislocations 
respectively, are such that pdis <Cpc^life^l- A 
moderate increase of e can cause a rapid increase of Pc up 
to its saturation value, while pdis may still remain very 
low. This situation can happen in the range of quench 
depths defined by {aScY < e < {o-^dis)^ (o is introduced 
in Ea. ((T^ 'l. provided that this range is sufficiently broad. 
Hence, one would expect the coarsening dynamics to be 
relatively insensitive to the value of e in that range, as ob- 
served numerically for 5|rc| < e < 35|rc| in Fig. 1111 Fur- 
ther increase of e eventually produces dislocation pinning, 
and a general slowing down of the system must occur: it 
is illustrated by the infiection of the effective exponents 
past 35|rc|. 

Let us assume next that the system is described by two 
characteristic length scales, Ldis and Lc, representing the 
linear extent of a grain along the y and a; directions, re- 
spectively. Following the arguments of [l^l, the dislo- 
cation density can be written as pdis ^ Ldis/ {LdisLc) = 
L^^. In the moderate quench regime (Sjrd < e < 35| rd), 
Lc ^ t^/'^. In this regime, the results on the energy 
and the structure factor fFig lll|l suggest that the other 
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length scale, Ldis, grows faster than t^/'^ (2 < z* < 3). 
Therefore, as time goes, defected regions tend to be 
more composed of dislocations than chevron boundaries 
{Ldis > Lc). Unfortunately, this finding disagrees with 
the experimental results on electroconvection, where the 
opposite behavior was found: a very slow growth law for 
Ldis (~ t^^^), as well as a similar law for a correlation 
length, was reported in j^. Hence, the experimental 
gra ins are elongated along the x direction at late stages 

V. CONCLUSIONS 

We have presented evidence that the coarsening of 
smectic patterns, as given by a potential anisotropic 
Swift-Hohenberg equation, is characterized by a t^^^ law 
close to onset. This law has not been observed exper- 
imentally so far in electroconvection of nematic liquid 
crystals, and may correspond to a regime difficult to 
reach. For larger quench depths, the phase ordering ki- 
netics is affected by pinning effects that strongly depend 



on grain boundary orientations. A particular regime is 
observed numerically for a fairly wide interval of moder- 
ate quenches: chevron boundaries get pinned, and grain 
growth still takes place via mobile arrays of dislocations. 
A similar behavior was observed in experiments ,.14„ „15| . 
In this regime, the dislocation density behaves as t~^^^, 
the same decay rate as found experimentally |15|. Our 
results suggest that the characteristic length of a disloca- 
tion array grows faster than that of a chevron boundary. 
This feature disagrees with the experimental findings, 
though, and may points out a limitation of the present 
model. Nonpotential effects, that have been neglected 
here, probably play an important role in this system. 
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